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The detection of fast variations of the TeV (10 12 eV) 7-ray flux, on time-scales 
of days, from the nearby radio galaxy M 87 is reported. These variations are 
~ 10 times faster than that observed in any other waveband and imply a very 
compact emission region with a dimension similar to the Schwarzschild radius 
of the central black hole. We thus can exclude several other sites and processes 
of the 7-ray production. The observations confirm that TeV 7-rays are emitted 
by extragalactic sources other than blazars, where jets are not relativistically 
beamed towards the observer. 

So far the only extragalactic objects known to emit 7-radiation up to energies of Tera elec- 
tron volts (1 TeV = 10 12 eV) are blazars. These are active galactic nuclei (AGN) with a plasma 
jet emanating from the vicinity of the black hole and pointing close to the observer's line of 
sight. Due to the bulk relativistic motion of the plasma in the jet the energy and luminosity of 
emitted photons are boosted by relativistic effects, making blazars detectable up to TeV ener- 
gies. 

The nearby radio galaxy M 87 is located in the Virgo cluster of galaxies at a distance of 
~ 16Mpc (z = 0.0043) and hosts a central black hole of (3.2 ± 0.9) x 10 9 solar masses 
©. The 2 kpc scale plasma jet © originating from the centre of M 87 is resolved at different 
wavelengths (radio, optical and X-rays). The observed inclination of the jet, at an angle of ~ 30° 
relative to the observer's line of sight (3 ), demonstrates that M 87 is not a blazar and hence would 
represent a new class of TeV 7-ray emitters. M 87 has also been suggested as an accelerator of 
the enigmatic ultra-high-energy (10 20 eV) cosmic rays (HE). Previously, weak evidence for E > 
730 GeV 7-ray emission from M 87 in 1998/1999 with a statistical significance of 4.1 standard 
deviations was reported by the High Energy Gamma Ray Astronomy (HEGRA) collaboration 
©. No significant emission above 400 GeV was observed by the Whipple collaboration © 
from 2000-2003. 

The observations reported here were performed with the High Energy Stereoscopic System 
(H.E.S.S.) located in Namibia. H.E.S.S. is an array of four imaging atmospheric-Cherenkov 
telescopes used for the measurement of cosmic 7-rays of energies between 100 GeV and several 
10 TeV, see © for more details. The observations of M 87 were performed between 2003 and 
2006 yielding a total of 89 hours of data after quality selection cuts. After calibration (|52T) , 
the H.E.S.S. standard analysis was applied to the data using hard event selection cuts (TSTb . 
More information about the standard analysis, as well as a more recent, alternative analysis 
technique (15*31) which gives consistent results, can be found in ((721) . 

An excess of 243 7-ray events is measured from the direction of M 87 in the whole data- 
set, corresponding to a statistical significance of 13 standard deviations, establishing M 87 as a 
TeV 7-ray source (Fig. 1). The position of the excess (Right Ascension a and Declination 5) 
was found to be a = 12 h 30 m 47.2 s ± 1.4 s , 5 = + 12°23'51" ± 19" (J2000.0). This is, within 
the quoted statistical error and the systematic pointing uncertainty of the H.E.S.S. telescopes 
(~ 20" in both the Right Ascension and Declination directions) compatible with the nominal 
(radio) position (0 of the nucleus of M 87 (a = 12 h 30 m 49.4 s , 5 = +12°23'28"). Considering 
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the angular resolution of H.E.S.S., the source is consistent with a point-like object with an 
upper limit for a Gaussian surface-brightness profile of 3 arcmin (99.9% confidence level). At 
the distance of M 87 (16 Mpc) this corresponds to a radial extension of 13.7kpc which can be 
compared with the large-scale structure of M 87 as seen at radio wavelengths (Fig. 1). A ~ 10 6 
times stronger constraint on the size of the TeV emission region is deduced from the observed 
short-term flux variability, as shown below. 

The differential energy spectra obtained for the 2004 and 2005 data sets (Fig. 2) are both 
well fit by a power-law function dN/dE oc E~ T . The spectrum measured in 2005 is found to 
be hard (r ~ 2.2) and reaches beyond 10 TeV with an average 7-ray flux of a factor of ~ 5 
higher than in 2004. 

The total 7-ray flux above 730 GeV (Fig. 3) for the individual years from 2003 to 2006 in- 
dicatates variability on a yearly basis (14\ corresponding to a statistical significance of 3.2 stan- 
dard deviations, being derived from a x 2 fit of a constant function. The variability is confirmed 
by a Kolmogorov test comparing the distribution of photon arrival times to the distribution 
of background arrival times yielding a statistical significance for burst-like (non-constant) be- 
haviour of the source of 4.5 standard deviations. Surprisingly variability on time-scales of days 
(flux doubling) was found in the high state data of 2005 (Fig. 3, upper panel) with a statisti- 
cal significance of more than 4 standard deviations. This is the fastest variability observed in 
any waveband from M 87 and strongly constrains the size of the emission region of the TeV 
7-radiation, which is further discussed below. No indications for short-term variability were 
found in the data of 2003, 2004 and 2006, which is not unexpected given the generally lower 
statistical significances of the 7-ray excesses in those years. 

These observational results (location, spectrum & variability) challenge most scenarios of 
very-high-energy 7-ray production in extragalactic sources. Although the luminosity (ps 3 x 
10 40 erg/s) of TeV 7-rays is quite modest and does not cause any problems with the global 
energy budget of the active galaxy M 87, several models can be dismissed. The upper limit 
on the angular size of ~ 3 arcmin (13.7kpc w 4.3 ■ 10 22 cm) centred on the M87 nucleus 
position already excludes the core of the Virgo cluster (15) and outer radio regions of M 87 
as TeV 7-ray emitting zones. Further, the observed variability on time-scales of At ~ 2 days 
requires a very compact emission region due to the light-crossing time. The characteristic size 
is limited to R < c ■ At ■ 5 ~ 5 x 10 15 5 cm « 5 x 5 R s , where 5 is the relativistic Doppler 
factor (17751) of the source of TeV radiation and R s ps 10 15 cm is the Schwarzschild radius of the 
M 87 supermassive black hole (see below). For any reasonable value of the Doppler factor (i.e. 
1 < 5 < 50, as used in the modelling of TeV 7-ray blazars), this implies a drastic constraint 
on the size of the TeV 7-ray source which immediately excludes several potential sites and 
hypotheses of 7-ray production. First of all this concerns the elliptical galaxy M 87 (1731) and the 
7-ray production due to dark matter annihilation (17). The most obvious candidate for efficient 
particle acceleration (18), namely the entire extended kiloparsec jet, is also excluded. Although 
compatible with the TeV source position, even the brightest knot in the jet (knot A) appears 
excluded with its typical size of the order of one arcsec (about 80 pc ps 2.5 ■ 10 20 cm) resolved 
in the X-ray range d7"9i 
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An interesting possibility would be the peculiar knot (HST-1) in the jet of M 87 (see support- 
ing online text and Fig. S2), a region of many violent events, with X-ray flares exceeding the 
luminosity of the core emission §S4\i and super-luminal blobs being detected downstream. Mod- 
elling the high-energy radiation properties of this region (by synchrotron and inverse-Compton 
scenarios), several authors favour sizes in the range of 0.1 to lpc (for moderate values of the 
Doppler factor ranging between 2 and 5) (|2ill22l|S?b . But, formally there is no robust lower limit 
on the size of HST-1, therefore we cannot exclude HST-1 as a source of TeV 7-rays. However, 
it would be hard to realize the short-term variability of the TeV 7-ray emission in relation to 
HST-1, at least within the framework of current models. While the size of the 7-ray production 
region does not exceed R < 5 x 10 15 5cm, the location of HST-1 along the jet at 0.85 arcsec 
from the nucleus, which corresponds to d ~ 65 pc, implies that the energy would be channelled 
from the central object into the 7-ray production region within an unrealistically small opening 
angle ~ R/d « 1.5 x 10~ 3 5 degree. 

The only remaining and promising possibility is to conclude that the site of TeV 7-ray 
production is the nucleus of M 87 itself (1231) . In contrast to the established TeV 7-ray blazars, 
the large scale jet of M 87 is seen at a relatively large jet angle (9 ~ 30°) which suggests a quite 
modest Doppler boosting of its radiation. Nevertheless, due to the proximity of M 87, both 
leptonic (24) and hadronic £5]|25) models predicted detectable TeV 7-ray emission. However, 
these scenarios typically produce a soft energy spectrum of TeV 7-rays, clearly in contrast to 
the hard spectrum measured by H.E.S.S. Leptonic models can be adapted in various ways to 
match the new results. Within synchrotron- self-Compton (SSC) scenarios (l26|) . one method 
is to consider the possibility of differential Doppler-boosting in the jet near the core region, 
a phenomenon clearly expected in the jet formation zone which extends over less than 0.1 pc 
from the nucleus (1271) . Emitting plasma blobs of small sizes with Doppler factors between 
5 and 30 and magnetic fields well below equipartition can account for the observed TeV 7- 
ray emission. An additional flux contribution from inverse-Compton scattering of background 
photons, coming from scattered disk emission or from dust, can further reduce the range of 
Doppler factors towards moderate values. 

The TeV 7-ray photons (independent of their production mechanism) might be absorbed by 
the pair absorption process 7 r p e y + 7jr — > e + e~ on the local infrared (IR) radiation field in 
the TeV 7-ray emission region. Since no signature for an absorption can be identified in the 
energy spectrum up to 10 TeV, one can derive an upper limit on the luminosity of the infrared 
radiation field at 0.1 eV (corresponding to a wavelength of approximately 10 micron, most 
relevant for absorption of 10 TeV 7-rays) to be L(O.leV) < 3.6 x 10 38 (i?/10 15 cm) erg/s, 
where R is the size of the TeV 7-ray emission region. Such a low central IR radiation luminosity 
supports the hypothesis of an advection-dominated accretion disk (i.e. an accretion disk with 
low radiative efficiency) in M 87 (28) and generally excludes a strong contribution of external 
inverse-Compton emission on IR light to the TeV 7-ray flux. 

If one accepts the hypothesis that protons can be accelerated as high as 10 20 eV in jets of ra- 
dio galaxies, then (hadronic) proton synchrotron models d5]|25l) can not be excluded considering 
the presented data. An alternative 7-ray production mechanism is curvature radiation of ultra- 
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high-energy protons in the immediate vicinity of the supermassive black hole. This novel mech- 
anism can simultaneously explain both the hard spectrum and fast variability of the observed 
TeV 7-ray emission. Rapidly rotating black holes embedded in externally supported magnetic 
fields can generate electric fields and accelerate protons to energies up to 10 20 eV (|29ll30ll3il) . 
Assuming that acceleration of protons takes place effectively within 3 Schwarzschild radii R s , 
and if the horizon threading magnetic field is not much below 10 4 G, one should expect 7-ray 
radiation due to proton curvature radiation extending to at least 10 TeV (the electron curvature 
radiation is less likely because of severe energy losses even in a tiny component of an irregu- 
lar magnetic field). No correlation with fluxes at other wavelengths is expected in this model. 
Although the size of the 7-ray production region, R ~ 3R S ~ 3 x 10 15 cm perfectly matches 
the observed variability scale, and the model allows extension of the 7-ray spectrum to 10 TeV 
without any significant correlation at other wavelength, the main problem of the model is the 
suggested magnetic field. It is orders of magnitude larger than the B-field expected from the 
accretion process, given the very low accretion rate as it follows from the bolometric luminosity 
of the core as well as the estimates of the power of the jet in M 87. 

In summary, the time-scale of the short-term variability of the TeV 7-rays is in the order of 
the light crossing time of the black hole (located at the center of M 87), which is a natural time- 
scale of the object. Therefore, the results reported here give clear evidence for the production 
of TeV 7-rays in the immediate vicinity of the black hole of M 87. 
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Fig. 1. Shown are the sky map as well as the position and extension limit of the TeV 7-ray emis- 
sion from M 87. (A) Smoothed TeV 7-ray excess map (color coded, 0.1° integration radius) as 
measured by H.E.S.S. The size (68% containment radius) of the H.E.S.S. point spread func- 
tion (PSF) is also indicated. The red circle indicates the intrinsic extension upper limit (99.9% 
confidence level) of 3 arcmin of the TeV 7-ray excess corresponding to 13.7kpc in M 87. The 
contour lines show the 90 cm radio emission (1321) . The white box marks the cut-out shown in the 
right panel. (B) The 90 cm radio data (1321 measured with the Very Large Array, together with 
the TeV position with statistical and 20" pointing uncertainty errors (white cross) and again the 
99.9% c.l. extension upper limit (red circle). Note that the size of the emission region deduced 
from the short-term variability is ~ 10 6 times smaller. The black cross marks the position and 
statistical error of the 7-ray source reported by HEGRA. The green ellipse indicates the host 
galaxy seen in the optical wavelengths with an extension of 8.3 x 6.6 arcmin in diameter. 
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Energy (TeV) 

Fig. 2. The differential energy spectrum of M 87 obtained from the 2004 and the 2005 data 
(using standard event selection cuts (1571) ) covering a range of ~ 400 GeV to ~ 10 TeV. Spectra 
for the 2003 and 2006 data sets could not be derived due to limited event statistics. Flux points 
with a statistical significance less than 1.5 standard deviations are given as upper limits (99.9% 
c.L). The corresponding fits of a power-law function dN/dE — I ■ (E/l TeV) _r are indicated 
as lines. The photon indices are T = 2.62 ± 0.35 (2004 data) and T = 2.22 ± 0.15 (2005 
data). Aside from the difference in the flux normalisation by a factor of ~ 5 (I = (2.43 ± 
0.75) x 10~ 13 cm" 2 s- 1 TeV" 1 in 2004 and J = (11.7 ± 1.6) x 10~ 13 cm- 2 s" 1 TeV" 1 in 2005) 
no variation in spectral shape is found within errors. The systematic error on the photon index 
and flux normalisation are estimated to be Ar = 0.1, and AI /I = 0.2, respectively. 
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Fig. 3. Gamma-ray flux above an energy of 730 GeV as a function of time. The given error 
bars correspond to statistical errors. (B) The average flux values for the years 2003 to 2006 as 
measured with H.E.S.S. together with a fit of a constant function (red line). The flux reported 
by HEGRA is also drawn (a systematic error must be taken into account when comparing re- 
sults from the two instruments). (A) The night-by-night fluxes for the four individual months 
(February to May) of the high-state measurements in 2005 with significant variability on (flux 
doubling) time-scales of ~ 2 days. The green points in (B) correspond to the 0.2 — 6 keV X-ray 
flux of the knot HST-1 (solid, ftS4\i ) and the nucleus (dashed, (TJ3l) ) as measured by Chandra; 
the lines are linear interpolations of the flux points. No unique correlation between the flux of 
X-rays and TeV 7-rays can be identified (the X-ray/TeV data were not taken simultaneously). 
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E-mail: matthias .beilicke @ desy.de, Olivier.Martineau-Huynh @ Ipnhep.in2p3 .fr 
H.E.S.S. observations and standard analysis 

When a primary 7-ray photon (or hadron) of TeV energies enters the earth's atmosphere, an 
extended air shower consisting of millions of secondary particles develops in the atmosphere. 
The Cherenkov light which is emitted in the air shower is recorded by the photomultiplier cam- 
era situated in the focal plane of each H.E.S.S. Cherenkov telescope. The recorded images 
of the shower are parameterized using a set of parameters (referred to as Hillas parameters) 
and are used to reconstruct for example the energy and direction of the primary particle. The 
geometrical width of an image is scaled using an expectation value for the corresponding obser- 
vation conditions (zenith angle, shower distance, etc.) and is used to distinguish between 7-ray 
and hadron induced showers. This procedure is referred to as the H.E.S.S. standard analysis 
and is described in more detail elsewhere (1571) . The simultaneous recording with up to four 
telescopes (stereoscopic observation) with this new generation experiment allows an improved 
(as compared to single telescope observations) event-by-event measurement of the direction 
(A9 ~ 0.1°) and energy (AE/E « 15%) as well as a superior background suppression of 
charged cosmic rays. Hard event selection cuts have been used in the analysis reported here 
providing an optimum background rejection for sources with hard energy spectra at some ex- 
pense of a lower event rate and an increased energy threshold of 640 GeV for the average zenith 
angle of 40° of the M 87 observations. The energy spectrum was derived using the standard cuts 
with a lower threshold of ~ 400 GeV for the same zenith angle. 

The H.E.S.S. observations of M87 were performed between 2003 and 2006 for a total of 
89 hours after data quality selection cuts (Tab. SI). The 2003 data (~ 25 h) were taken with 
only two operational telescopes during the construction phase, while the measurements of the 
following years have been performed with the full four telescope array which has a twice better 
sensitivity as compared to the two-telescope setup. Therefore, the energy spectra and the sky 
position and extension limit of the TeV 7-ray excess were derived from the 2004-2006 data 
only. Generating an energy spectrum from the 2006 data was not possible due to the very 
limited event statistics. 

Application of an alternative analysis method 

The results have been cross-checked with a recently developed alternative analysis method. Be- 
side an independent calibration of the H.E.S.S. raw data (15*21) , this technique is based on the 
combination of a parameterisation of the shower images using the moment method of Hillas 
(similar to the one described above), and the technique refered to as model analysis (lS3l) . This 
latter method implies a pixel-by-pixel comparison of the shower images (recorded with the 
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H.E.S.S. photomultiplier cameras) with a template generated by a semi-analytical shower de- 
velopment model. The 7-ray primary energy, direction and impact position are obtained by 
maximising a log-likelihood function associated with this comparison. The model analysis uses 
all available pixels in the camera, without the requirement of an image cleaning. Since the 
background rejection is based on independent variables in these two methods, their combina- 
tion in this alternative analysis improves the background rejection. The combined analysis used 
here therefore yields more significant results than the Hillas parameter based standard analysis 
alone. 

The light curve of the 2005 flux high state of M 87, in which the short-term variability is 
identified, obtained with this alternative analysis method is shown in Fig. S 1 . The corresponding 
significance of the short-term variability is above 6 standard deviations, clearly confirming the 
results obtained from the standard analysis described in the main text. The energy spectra (T = 
2.61 ±0.24, J = (3.60 ± 0.57) x 1(T 13 cm~ 2 s^TeV^ 1 for the 2004 data and T = 2.20 ±0.09, 
I Q = (13.9 ±1.2) x 10~ 13 cm" 2 s" 1 TeV 1 for the 2005 data) as well as the sky position of 

a = 12 h 30 m 50.0 s ± 1.3 s and 5 = ±12°23'53" ± 19" (J2000.0) and extension limit of 3 arcmin 
of the TeV 7-ray excess are also compatible with the results of the standard analysis. 

The peculiar knot HST-1 

HST-1 is the innermost resolved knot in the jet of M 87 (Fig. S2). Up to now, no extension or 
substructure of HST-1 could be resolved in the X-ray range. This region is the site of many vi- 
olent events, with X-ray flares exceeding the luminosity of the core emission and super-luminal 
blobs being detected downstream. The radio, UV and X-ray fluxes ($4) of HST- 1 have increased 
by more than a factor of 50 over the period from 2000 to 2005 (Fig. 3). In addition, isolated 
flares with variability on time- scales of about 1 month have been detected both in optical (TS3T) 
and X-ray (IS?!) energy bands. 
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Fig. SI. Gamma-ray flux above an energy of 730 GeV as a function of time for the 2005 
high state (February to May) from the two analysis methods in comparison. The given error 
bars correspond to statistical errors. (A) The results of the standard analysis shown in Fig. 3, 
whereas (B) shows the fluxes obtained using the model analysis (lS3l) , including an independent 
calibration chain. Due to its higher 7-ray acceptance, about 60% of the 7-ray events obtained 
from the model analysis are not contained in the corresponding standard analysis in which the 
hard event selection cuts were used, making the two data sets (partly) statistically independent. 
Both methods exhibit short-time variability with comparable time-scales. 
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Fig. S2. The central 2 kpc jet of M 87 which is resolved at radio (A), optical (B) and X-ray (C) 
energies. The nucleus (N), the closest knot HST-1 as well as other knots in the jet are indicated. 
The image is adopted from (S6). 
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Tab. SI. Summary of H.E.S.S. observations of M 87. Shown is the time period for each year 
in which observations were performed, the observation time T Q ^ S and Tij ve (corrected for the 
detector dead time) as well as the number of telescopes in the system. The 2003 data (marked 
with a *) were taken with only two telescopes in a single operation mode (no stereoscopic 
trigger) during the construction phase. 



Year 


time period 


Tobs M 


Tlive M 


^tel 


2003* 


24. Apr. - 28. June 


25.1 


19.2 


2 


2004 


17. Feb. - 23. May 


35.2 


31.4 


4 


2005 


12. Feb. - 15. May 


23.0 


21.1 


4 


2006 


25. Apr. - 18. May 


5.3 


4.8 


4 
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